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Abstract 

We derive high-temperature series expansions for the free en- 
ergy and susceptibility of the two-dimensional random-bond Ising 
model with a symmetric bimodal distribution of two positive coupling 
strengths J\ and J2 and study the influence of the quenched, random 
bond-disorder on the critical behavior of the model. By analysing the 
series expansions over a wide range of coupling ratios J2/J1, covering 
the crossover from weak to strong disorder, we obtain for the suscepti- 
bility with two different methods compelling evidence for a singularity 
of the form x ~ t~ 7 ^ 4: \ lni| 7//8 , as predicted theoretically by Shalaev, 
Shankar, and Ludwig. For the specific heat our results are less con- 
vincing, but still compatible with the theoretically predicted log-log 
singularity. 



1 Introduction 



One of the most-studied variants of the two-dimensional Ising model is the 
case of random bonds. While realizations of Ising models that include ran- 
domness come much closer to approximating reality, they are very much 
harder to study at any level. Even in two dimensions exact results for random 
cases (especially for quenched randomness, which is the realistic situation in 
many experiments) are few and far-between. In fact, the two-dimensional 
Ising case is especially difficult because of the marginality of the Harris crite- 
rion U for this model. This criterion states that quenched randomness is a 
relevant (irrelevant) perturbation when the critical exponent a of the specific 
heat of the pure system is positive (negative) and therefore when a = as 
in the two-dimensional Ising model the situation is marginal. 

Numerous theoretical investigations [2-7] as well as numerical Monte 
Carlo simulations [8-15] and transfer- matrix studies J16, [17| have addressed 
the question of whether the critical exponents for the two-dimensional Ising 
model with quenched, random bond disorder differ from those of the pure 
model. While a "majority" consensus had probably been achieved in favour 
of no change, apart from logarithmic corrections [3-6] no unambiguous nu- 
merical study that confirmed the quantitative predictions of either of the 
theoretical approaches had been made prior to our recent study of the sus- 
ceptibility with high-temperature series expansions. In a brief note fl8| we 
announced the confirmation of the theoretical majority consensus value of 
the exponent of the logarithmic correction. This quantitative determina- 
tion of the value of the correction exponent in excellent agreement with the 
predicted value, using a completely different numerical approach that in no 
way depends on random numbers, provided additional and incontrovertible 
support to the previous consensus. In the present paper we present the coef- 
ficients of the susceptibility series that we analysed []TB[ together with some 
remarks on their derivation, the details of our analysis, and some additional 
results for the specific-heat series. We note that a recent simulation of the 
site-diluted model [pl| also confirms the log- log prediction of [3-6] . 

In the next section, we define the model and the quantities that are 
studied, and in Sec. 3 the theoretical predictions are briefly recalled. The 
series generation is described in Sec. 4, and in Sec. 5 we describe the analysis 
techniques used. Sec. 6 then presents the results, where details of our analyses 
for the susceptibility give compelling evidence for a singularity of the form 



1 



predicted by Shalaev, Shankar, and Ludwig [3-6]. In Sec. 7 we close with a 
summary of our conclusions and a few final comments. 



2 Model 

The Hamiltonian of the random-bond Ising model is given by 

H = -^J ij a i a j , (1) 

where the spins er» = ±1 are located at the sites of a square lattice, the symbol 
(ij) denotes nearest-neighbor interactions, and the coupling constants J^ 
are quenched, random variables. As in most previous studies we consider a 
bimodal distribution, 

P{J tj ) = xSiJij - JO + (1 - x)S(J ij - JO, (2) 

of two ferromagnetic couplings Ji, J 2 > 0. We furthermore specialize to a 
symmetric distribution with x = 1/2, since in this case the exact critical 
temperature T c can be computed for any positive value of J\ and J2 from the 
(transcendental) self-duality relation (ks = Boltzmann constant) |20| 



(exp(2J 1 /A; B T c ) - 1) (exp(2J 2 /£; B T c ) - 1) = 2. (3) 

In both, computer simulations and high-temperature series expansion stud- 
ies, this exact information simplifies the analysis of the critical behavior 
considerably. 

The free energy per site is given by, 



Pf=- Km ^ 
V—+00 y 



In II E exp(-/W) 



ii=l CTi=±l, 



(4) 



where (3 = l/fc^T is the inverse temperature in natural units and the bracket 
[. . .] av denotes the average over the quenched, random disorder, [. . .] av = 
I dJi^j (. . .)P(Jij). The internal energy and specific heat per site follow 
by differentiation with respect to f3, 

e = dPf/dp, C/k B = -(3 2 d 2 (3f/d(3 2 . (5) 
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In this paper we shall mainly focus on the magnetic susceptibility per site \ 
which in the high-temperature phase and zero external field is defined as the 
V — > oo limit of 

(6) 

where (. . .)t denotes the usual thermal average with respect to exp(— PTC). 




3 Theoretical predictions 

Let us briefly recall two contradicting theoretical predictions for the critical 
behavior of the model (H), @. The first is based on renormalization-group 
techniques developed by Dotsenko and Dotsenko (DD) 0. For the specific 
heat DD find close to the transition point a double logarithmic behavior, 

C(t) cc ln(ln(l/|t|)), (7) 

where t = (T—T c ) /T c denotes the reduced temperature. For the susceptibility 
they obtain in the high-temperature phase (t > 0) 



X oc t 2 exp 



—a ( I n hi ( — 



(8) 



The second approach by Shalaev, Shankar, and Ludwig (SSL) [3-6] makes use 
of bosonisation techniques and the method of conformal invariance. While 
the prediction (|7|) for the specific heat can be reproduced (which, however, 
is not undisputed 0), SSL derive quite a different behavior for the suscepti- 
bility, 

XOct- 7/4 |hit| 7/8 . (9) 

This is the same leading singularity as in the pure case (J\ = J2), but modi- 
fied by a multiplicative logarithmic correction. 

High-precision Monte Carlo simulations and transfer- matrix studies [8- 
16] favor the latter form, but due to well-known inherent limitations of this 
method it has been impossible to confirm the value of the exponent of the 
multiplicative logarithmic correction in (|9]) quantitatively. Similar problems 
have been reported in simulation studies of other models exhibiting multi- 
plicative logarithmic corrections such as, e.g., the two-dimensional 4-state 
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Potts |21| and XY |2^] model. We found it therefore worthwhile to investi- 
gate this problem yet again with an independent method. In the following 
we report high-temperature series expansions for the free energy and suscep- 
tibility and enquire whether series analyses can yield a more stringent test 
of the theoretical predictions. 

4 Series generation 

For the generation of the high-temperature series expansions of the free en- 
ergy (HD, and hence the internal energy and specific heat, as well as the 
infinite- volume limit of the susceptibility (^) we made use of a program pack- 
age developed at Mainz originally for the g-state Potts spin-glass problem 
[23-27] In this application the spin-spin interaction is generalized from a^aj 
to 6 ai;aj with <7j being an integer between 1 and q, and the coupling constants 
Jij can take the values ± J with equal probability. Since here the coupling 
constants also can take negative values, frustration effects play an important 
role and the physical properties of spin glasses |28j are completely different 
than those of the random-bond system. Technically, however, precisely the 
same enumeration scheme for the high-temperature graphs can be employed 
in both cases. The only difference is in the last step where the quenched 
averages over the are performed. The details of the employed star-graph 
expansion technique and our specific implementation are described elsewhere 
[23-25, 27, 29] Here we only note that slight modifications of this program 
package enabled us to generate the high-temperature series expansion for ftf 
and x U P t° t ne 11th order in k = 2(3 J\ for 

• hypercubic lattices of arbitrary dimension d, 

• arbitrary number of Potts states q, 

• arbitrary probability x in the bimodal distribution, and 

• arbitrary ratios R = J2/J1, characterizing the strength of the disorder. 

In this paper we shall concentrate on the two-dimensional (d = 2) random- 
bond Ising model (q = 2) for a symmetric (x = 1/2) bimodal distribution 
of two positive coupling strengths J\ and J 2- The series coefficients of the 
free energy and susceptibility expansions for various coupling-strength ratios 
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R = Jij J\ are given in Tables [L] and |2|. Of course, in principle it would be 
also straightforward to adapt the present program package to more general 
probability distributions P(Jij). 



5 Series analysis techniques 

In the literature many different series analysis techniques have been discussed 
which, depending on the type of critical singularity at hand, all have their 
own merits and drawbacks In the course of this work we have tested 



quite a few of them [ 29| . Here, however, we will confine ourselves to only 
those which turned out to be the most useful for our specific problem at 
hand. 

To simplify the notation we denote a thermodynamic function generically 
by F(z) and assume that its Taylor expansion around the origin is known up 
to the iV-th order, 

N 



F(z) = Y.a n z n + .... (10) 

n=0 

If the singularity of F(z) at the critical point z c is of the simple form (z < z c ) 

F(z)~A(l-z/z c )-\ (11) 

with A being a constant, then the ratios of consecutive coefficients approach 
for large n the limiting behavior 



r n = 



dn-l 



A 



n 



1 

Z, 



(12) 



From the offset (l/z c ) and slope ((A — l)/z c ) of this sequence as a function of 
1/n both the critical point z c and the critical exponent A can be determined. 
This is the basis of the so-called ratio method [31]. If the critical point z c 



is known from other sources (in our case exactly from self-duality), then one 
may consider biased extrapolants for the critical exponent, 



A, 



n 



(13) 



which simply follow by rearranging eq. fll2). In the following this method 
will be denoted as "Biased Ratio I" . 
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If the singularity of F(z) contains a multiplicative logarithmic correction 
(as, e.g., in the SSL prediction for x), 

F(z) ~ A(l - z/z c )~ x \ ln(l - z/z c )\ p , (14) 

then one forms the ratios r n as before, but considers in addition the auxiliary 
function |32| 



Z 



-A 



(ln[l/(l - *)])*• = £ &„*" + ..., (15) 



n=0 



and computes the ratios r* = b n /b n -\. Let us first assume that the critical 
exponent A of the leading term is known. Then it can be shown that the 
sequence R n = r n jr* n approaches l/z c with zero slope in the limit n — > oo, 
if and only if p* = p. This determines p, if also z c is known. If A is not 
known, then one may vary both exponents until above relation is satisfied. 
In the following we refer to this special ratio method as "Ln-Ratio" . 

Another method [[3^, |33J suitable for a singularity of the form ( jTjj ) is 
based on Pade approximants | 55| . Here one generates the series expansion 
for the auxiliary function 

G(z) = -{z c - z) ln(z c - z)(^ - (16) 

which can easily be shown to satisfy 

lim G{z) = p. (17) 

If z c is known, the value of G(z) at z = z c can be obtained by forming Pade 
approximants, 



C(A „ rr /M] _ Pl{z) _ p + Pl Z+p 2 Z 2 + ...+p L Z L 

G{Z) ~ [L/M] = 0^) = l + qi z + q2Z > + ... + qMZ M> (18) 

with L + M < N — 1. Note that one order of the initial series is lost due to 
the differentiation in ([TBI ). 

With a small modification this method can also be applied to a purely 
logarithmic singularity of the form 

F(z) ~ A|ln(l - z/z c )\ p . (19) 
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In this case one defines the auxiliary function 



G(z) 



(z c - z) ln(z c - z) 



F\z) 



(20) 



F{zY 



which again satisfies 



lim G{z) = p. 



(21) 



The two analysis methods based on Pade approximants will be called "Ln- 
Pade". 

6 Results 



Susceptibility: In a first step we investigated whether our series expan- 
sions for the susceptibility are consistent with a pure power-law behavior 
according to the DD prediction (H) (ignoring the exponentially small multi- 
plicative correction term). Assuming thus the behavior \ oc t -7 and using 
the method "Biased Ratio I" we obtained the critical exponents 7 shown in 
Fig. [I] as a function of J^j ' J\- Here and in the following the error bars are 
estimated by varying the length of the series and/or the type of Pade approx- 
imants used. Starting with 7 = 1.738 ± 0.014 for the pure case (J2/ J\ = 1), 
being consistent with the exact value of 7 = 7/4, we observe a steady increase 
to 7 = 2.37 ± 0.11 for the strongest investigated disorder (J2/J1 = 10). We 
will argue below that the apparent crossover from weak to strong disorder is 
due to the finite length of our series expansion which naturally has a much 
more dramatic influence for weak disorder. At any rate, for strong disorder 
the DD prediction of 7 = 2 is clearly outside the error margins of the series 
analysis estimates. 

So far no multiplicative logarithmic corrections were taken into account. 
If the SSL prediction (|) was correct we would, therefore, expect to observe 
"effective" critical exponents which according to 



should indeed be larger than 7/4. The results in Fig. [T] could thus be well con- 
sistent with a critical exponent of 7 = 7/4 in the presence of a multiplicative 
logarithmic correction. 



X oct- 7 / 4 |lnt| 7 / 8 = t 



(7/4) [1+1 



1 ln(| lnt|) 

2 ln(l/t) 



(22) 



7 



This possibility suggested a more careful analysis based on the qualitative 
form of the SSL prediction @. Our series are too short to employ a general 
ansatz with both exponents as free parameters. We rather fixed the exponent 
7 = 7/4 of the leading term to the (predicted) pure Ising model value and 
enquired if our series expansions are compatible with the ansatz 



X oc 



t' 7/4 \\nt\ p , (23) 



and p — 7/8. Employing the two special methods for this type of singularity 
described in Sec. |5] we obtained well converging results. The resulting esti- 
mates for the exponent p are shown in Fig. [| We see that the two methods 
yield consistent results which start in the pure case (J2/ J\ = 1) around p = 0, 
as they should do. With increasing disorder the estimates exhibit again an 
apparent crossover, until around J2/ J\ = 5 — 8 they settle at a plateau value 
in very good agreement with the theoretical prediction of p = 7/8. This is 
the main result of our series analysis. We claim that compared with previous 
methods this is thus far the clearest quantitative confirmation of the SSL 
prediction (|9|). 

As before we attribute the apparent crossover for intermediate strength 
of the disorder to the shortness of our series expansions, i.e., we interpret the 
crossover as an unavoidable artifact of high-temperature series expansion 
analyses and not as an indication that the exponent p really is a function of 
the disorder strength. We thus take the view that already a small amount 
of disorder drives the system into a new universality class different from 
the pure case which, however, only becomes visible in the very vicinity of 
the transition point T c (or t = 0). This in turn translates into the need of 
extremely long series expansions in order to be detectable. 

To justify this claim we have investigated a model function simulating 
the "true" susceptibility (g > 0), where g is a constant that depends on 
the strength of the disorder, 



Xmodcl t ^ 



l + ^ln(l/i) 

7T 



7/8 



(24) 



with t = (T — T c )/T, which for any g 7^ reproduces the SSL form (||) in 
the limit T — > T c (t = t — t 2 + t 3 + . . . — > 0). Notice the discontinuity in 
the asymptotic behavior at <?o = 0. For any go ^ the asymptotic region is 
reached when ln(l/t) is much larger than 7r/4g Q , i.e., for t exp(— 7r/4g ) • 



8 



Since go = corresponds to the pure case it is intuitively clear that the 
parameter go is an increasing function of the degree of disorder. For weak 
disorder this implies that go is very small and therefore, due to the exponen- 
tial dependence, that the asymptotic region in t is extremely narrow. 

Strictly speaking the model function fl24[) should resemble the "true" sus- 
ceptibility only for weak disorder, but it is commonly believed that it is a 
reasonable qualitative approximation also for strong disorder. To relate the 
parameter go at least heuristically to the ratio J2/J1 we used the weak dis- 
order result go = C2Q 1 /(l + ab) 2 , where C2 = 1 — x (with x as defined in 
eq. (0)) is assumed to be small, C2 < 1, i.e., the analytic calculation as- 
sumes that there are only few J2-bonds in a background of Ji-bonds. The 
parameters a and b are given by a = (v' c — ) /v^ and b = v^'/2y/2, where 
= tanh(/?(°)Ji) = V2 - 1 and v' c = tanh($°> J 2 ), with (3^ denoting the 
inverse critical temperature of the pure system with all J^- = J\. Of course, 
employing this formula to the present case with C2 = l/2 = xisa bold 
step which even creates an ambiguity since the exact symmetry J\ <-» J2 for 
x = 1/2 is violated. For weak disorder (J 2 /Ji ~ 1), however, the inconsis- 
tency turns out to be very mild. For J 2 / J\ = 1.2 we obtain g = 0.013700 . . ., 
and for J2/J1 = 1/1.2 we find a slightly smaller value of g' = 0.011958. . .. 
This shows that for weak disorder {Jij J\ = 1.2, go ~ 0.013) the asymptotic 
region is bounded by t <C exp(— vr/4g ) ~ exp(— 1/0.017) ~ 10 -26 , and thus 
explains why it is so difficult to observe the asymptotic critical behavior in 
the weak-disorder limit. For J2/J1 = 1.5 and 1/1.5 the corresponding num- 
bers are go = 0.070700 . . . and g' = 0.052889 . . ., leading to a bound of the 
order of t 10~ 5 . 

Using a symbolic computer program it is straightforward to generate the 
high-temperature series expansion of the model function (p4|) to any desired 
order. Applying the same analysis techniques as used for the "true" suscep- 
tibility series we obtained the results shown in Fig. [3[ If we truncate the 
model series at low order we observe qualitatively the same crossover effect 
as for the "true" series. Here we are sure, however, that this must be a pure 
artifact of the truncation of the model series at a finite order. We also see 
that the approach of the asymptotic limit of p = 7/8 as a function of the 
degree of disorder is faster if we consider a longer series (21 terms). It is, 
however, somewhat discouraging (even though understandable in view of the 
exponential dependence of the critical regime on go) that at a fixed go the 
convergence of the series with increasing order is quite slow. For example, at 
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4fi'o/ 7r = 1 we obtained p = 0.7056 with the Pade approximant [4/4], 0.7178 
([5/5]), 0.7474 ([10/10]), 0.7682 ([20/20]), 0.7777 ([30/30]), 0.7834 ([40/40]), 
0.7875 ([50/50]), and 0.7905 ([60/60]). The convergence behavior for this ex- 
ample and other small values of the parameter g can be visually inspected 
in Fig. g. 



Specific heat: Series analyses for the specific heat are usually more diffi- 
cult than for the susceptibility. This is especially pronounced for the Ising 
model on loose-packed lattices where all odd powers of (3 vanish because of 
symmetry. Consequently our specific-heat series consists only of four non- 
trivial terms (see Table [l]). We nevertheless tried an analysis with the ansatz 

Coc|lnt| 9 , (25) 

using the method "Ln-Pade" . The exponent q is an effective exponent whose 
value may, or may not be constant. 

The resulting dependence of the exponent q on the ratio J^j ' J\ is shown 
in Fig. |5|. While the quantitative agreement with the exactly known pure 
case is certainly not convincing, we do see at least a qualitative trend to 
smaller values of q with increasing strength of the disorder (increasing ratios 
J2/J1), i-e., the singularity of the specific heat becomes apparently weaker 
for stronger disorder. This may be taken as an indication that the true 
singularity is of the log-log type (0) as predicted by both, DD and SSL. A 
recent numerical study |L7 for J2/J1 = 4 using transfer-matrix methods also 



observed a behavior in between log and log-log type. These findings are 
in contradiction to the claim |36| for a slightly different disordered system 
(quenched, random site-dilution) that the specific heat stays finite at T c , as 
theoretically suggested in Ref. (see also Ref. [|3~7|] ). 



Again we have tried to justify our interpretation by considering a model 
function, 

1 



Cmodel — — hi 
90 



IT 



(26) 



By applying precisely the same type of analysis to the series expansion of the 
model specific-heat we obtained the results displayed in Fig. which show 
qualitatively the same trend of decreasing q as a function of J2/J1 as the 
data in Fig. H. 
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7 Discussion 



The main results of our high-temperature series analysis are shown in Fig. |2] 
which provide at least for strong disorder (large J2/J1) compelling evidence 
that the singularity of the susceptibility is properly described by x oc t~ 7 / 4 | lnt\ p , 
with p = 7/8 = 0.875, as theoretically predicted by SSL [3-6]. The analy- 
sis of the model susceptibility (0) in Figs. |3] and [| clearly shows that the 
apparent variation of p with the strength of disorder is an artifact caused 
by the truncation of the series expansions at a finite order. We, therefore, 
emphasize that the apparent crossover from weak to strong disorder does not 
imply that the universality class of the random-bond Ising model changes 
continuously with the strength of disorder. 

Let us finally make a few comments on previous Monte Carlo simula- 
tions of this model on large but finite square lattices. With the finite-size 
scaling analysis of Refs. [8-12, 16] it is conceptually impossible to detect the 
multiplicative logarithmic correction of the SSL prediction ([]). The reason 
is that the SSL theory also predicts a logarithmic correction for the scaling 
behavior of the correlation length, £ oc t" 1 ] hit I 1 / 2 . In the finite-size scaling 
behavior the two logarithms thus cancel and one ends up with a pure power- 
law, x °c = L 7 / 4 , where L is the linear lattice size. Thus only the SSL 
prediction for 7/1/ can be tested in finite-size scaling analyses. Wang et al. 
[|, obtained for J 2 j J x = 4 and 10 an estimate of y/v = 1.7507 ± 0.0014, 
and also the results of Reis et al. Jl6| at 3%j J\ = 2, 4, and 10 are consistent 
with 7/V = 1.75. Among the two alternatives, the theories of DD and SSL, 
these estimates thus provide evidence in favor of SSL. Notice, however, that 
a numerical estimate of 7/V ~ 1.75 would also be expected for the pure two- 
dimensional Ising model. For the specific heat the situation is conceptually 
clearer. Here the theoretically expected scaling behavior (0), as predicted 
by both, DD and SSL, translates into a double-logarithmic finite-size scaling 
behavior, C = Cq + C% ln(l + 61nL), which is different from that of the pure 
case where C = Cq + InL. In the numerical work of Wang et al. || [II]], 
employing lattice sizes up to L = 600, this difference in the asymptotic be- 
havior is clearly observed for J2/J1 = 10, while for J2/J1 = 4 the behavior is 
in between log and log-log type, similar to the findings reported in a recent 
transfer-matrix study [17] for the same coupling-constant ratio. For the spe- 
cific heat these latest finite-size scaling analyses are thus about as conclusive 
as our series analyses in Fig. |^. 
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Another set of numerical data that can discriminate between the pre- 
dictions of DD and SSL comes from direct simulations of the temperature 
dependence of the magnetization m and of the susceptibility x f° r J2/ Ji = 4 
|fL0| , [L3|]. Assuming in the analysis a pure power law with an effective expo- 
nent (i.e., ignoring the logarithmic correction), one observes an overshooting 
of the effective exponent to values larger than the prediction by SSL. As 
discussed above (recall eq. (|22|)) this may be taken as an indication of a mul- 
tiplicative logarithmic correction term. For example, Talapov and Shchur 
[[yl obtained for J%j J\ = 4 from least-squares fits to x °t t~ lcti an effective 
exponent of 7 e g ~ 7/4 + 0.135 = 1.885. This value is quite close to our series 
estimate of 7 e ff = 2.019 ± 0.024 for J2/J1 = 4, if the pure power-law ansatz 
is used (cp. Fig. |I]). Wang et al. |TtJ furthermore confirmed that their data 
is compatible with the SSL ansatz, x(t) = x t _7//4 (l + at)[l + £> ln( 1/ t)] 7//s , 
supplemented by a correction to scaling term (1 + at) (and similarly for m; 
for a recent confirmation, see Ref. [|HJ). In these fits both exponents are 
kept fixed at their predicted values, and xo, a , and b are free parameters. 
In contrast to our series analysis, however, no quantitative estimates of the 
exponent of the logarithmic correction have been reported in Ref. . While 
the simulation results certainly indicate that among the two conflicting theo- 
ries of DD and SSL, the SSL prediction is more likely to be correct, it is still 
fair to conclude that also this set of simulations has not yet unambiguously 
identified the multiplicative logarithmic correction term. 

Monte Carlo simulations of systems with quenched, random disorder re- 
quire an enormous amount of computing time because many realisations have 
to be simulated for the quenched average. For this reason it is hardly possible 
to scan a whole parameter range. Using high-temperature series expansions, 
on the other hand, one can obtain closed expressions in several parameters 
(such as the dimension d, x, J2/ J\, ■ ■ ■ ) up to a certain order in the inverse 
temperature f3 = l/fc^T. Here the infinite- volume limit is always implied 
and the quenched, random disorder can be treated exactly. By analysing the 
resulting series, the critical behavior of the random-bond system can hence 
in principle be monitored as a continuous function of several parameters. 
This is a big advantage over Monte Carlo simulations which usually can only 
yield a rather small parameter range in one set of simulations. The caveat of 
the series-expansion approach is that the available series expansions for the 
random-bond Ising model are still relatively short (at any rate much shorter 
than for pure systems). This introduces systematic errors of the resulting es- 
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timates for critical exponents which are difficult to control. The obvious way 
out is trying to extend the series expansions as far as possible. This, however, 
would be extremely cumbersome since the number of algebraic manipulations 
necessary to calculate the series coefficients blows up dramatically with the 
order of the series (usually at least exponentially) and, therefore, has to be 
left for future work. 
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Table 1: Coefficients a n of the high-temperature series expansion of the free 
energy per site, —f3f = In [2 cosh(/9Ji) cosh(/?J 2 )] + J2n a nk n , with k = 2/3 Ji. 



n 


J2/J1 = 1 


J2/J1 = 2 


J2/J1 = 4 


J2/J1 = 10 


4 


1 

16 


81 
256 


625 
256 


14641 

256 


6 


1 

96 


81 
2048 


18125 
6144 


5343965 
6144 


8 


17 
2560 


33671 
327680 


26293 
65536 


990135929 
327680 


10 


1907 


3437297 


1057390637 


16514750542133 


483840 


27525120 


49545216 


49545216 



Table 2: Coefficients b n of the high-temperature series expansion of the 
susceptibility per site, x — 1 + J2 n b n k n , with k = 2(3 J x . 



n 


J2/J1 = 1 


J2/J1 = 2 


J2/J1 = 4 


J2/J1 = 10 


1 


2 


3 


5 


11 


2 


3 


27 
4 


75 
4 


363 
4 


3 


13 

3 


231 
16 


3115 
48 


31933 
48 


4 


23 
4 


1809 
64 


13025 
64 


277937 
64 


5 


451 
60 


69337 
1280 


471185 
768 


101248147 
3840 


6 


191 
20 


515871 
5120 


1823875 
1024 


768499919 
5120 


7 


30283 
2520 


79576207 
430080 


1302083479 
258048 


1034056024661 
1290240 


8 


100003 
6720 


191638233 
573440 


4823704415 
344064 


7079050432267 
1720320 


9 


3318601 


587805509 


203928262469 


3850544162365417 


181440 


983040 


5308416 


185794560 


10 


3369629 
151200 


48645511629 
45875200 


5160699783175 
49545216 


126985060534491247 
1238630400 


11 


269543489 
9979200 


101837138460677 
54499737600 


9157142004160957 
32699842560 


1069481408075459203 
2123366400 
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Figure 1: Analysis of the susceptibility series assuming a singularity of the 
form x oc t" 7 7 using the method "Biased Ratio I". 
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Figure 2: Analysis of the susceptibility series assuming a singularity of the 
form x ^ t~ 7 ^ 4 \lnt\ p , using Pade approximants and the ratio method (see 
text). The horizontal line at p = 7/8 = 0.875 is the theoretical prediction of 
SSL. 
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Figure 3: Analysis of the series expansion of the model function using 
the ansatz Xmodci oc t~ 7 / 4 \ \nt\ p . The legend indicates the three different Pade 
approximants shown. The horizontal line shows the exact value p = 7/8 = 
0.875. 
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Figure 4: Convergence behavior of the model series for the susceptibility 
with increasing order at fixed parameter g . The symbols [L/M] denote the 
various Pade approximants used. 



18 



1.5 



0.5 



0.0 





<> 

- ^ $ 

<> § § § 












o Ln-Pade 


* 






* 



12 3 4 



J, /J, 

Figure 5: Analysis of the specific-heat series using the Ln-Pade method. 
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Figure 6: Analysis of the series expansion of the model specific-heat using the 
Ln-Pade method. 
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